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LOCALIZATION FOR THE RANDOM DISPLACEMENT MODEL 

FREDERIC KLOPP\ MICHAEL LOSS^ SHU NAKAMURA^, AND GUNTER STOLZ" 

Abstract 

We prove spectral and dynamical localization for the multi-dimensional random displace- 
O ■ ment model near the bottom of its spectrum by showing that the approach through multiscale 

analysis is applicable. In particular, we show that a previously known Lifshitz tail bound can 
be extended to our setting and prove a new Wegner estimate. A key tool is given by a quan- 
titative form of a property of a related single-site Neumann problem which can be described 
as "bubbles tend to the corners". 

1. Introduction 
P_il We consider the random displacement model (RDM), a random Schrodinger operator 

i : H^ = -A + V^ (1) 

-4— > . 

3 ' in L^(M'^), d > 1, where the random potential has the form 

V^{x) = J2(li^-^-'^i)- (2) 

J^ . This models a random perturbation of the periodic potential Yli q{^ ~ O^ where the single- 

00 \ site terms sit at exact lattices sites i G Z"'. The parameter u = {uJi)i^^d describes a configura- 

"^ I tion of random displacement vectors cui E M.'^. Before entering a more thorough discussion of 

background and assumptions, let us state the main result of our work in a simple non-trivial 
Q . special case: 



O 



^ 



Theorem 1.1. Suppose that d > 2 and that q G C^(R'^) is real-valued, sign-definite, rotation 
symmetric and suppq C {x : |a;| < r} for some r < 1/4. Also assume that {oJi)i(zzd are i.i.d. 
W^-valued random variables, uniformly distributed on [—dmaxidmaxY where dmax = 1/2 — ^■ 
Then H^ is spectrally and dynamically localized at energies near the bottom of its almost 



ci I sure spectrum. 



The exact meaning of the latter will be recalled below. 

The RDM represents a natural way to model a solid with structural disorder. It can be 
considered as intermediate between the Anderson model, which has no structural disorder and 
randomness instead appears in the form of coupling constants at the single site terms, and the 
Poisson model, where the structure of the medium is entirely dissolved by placing single-site 
scatterers at the points of a Poisson process. This point of view has recently been supported 
by an investigation of the integrated density of states of the RDM in [9] and [10]. 
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It is physically expected that multi-dimensional random Schrodinger operators such as the 
three models mentioned above should exhibit localization at energies near the bottom of the 
spectrum, while extended states should exist at high energy. For mathematicians the latter 
remains an open problem, while localization at low energy for the Anderson model (both in 
the continuum and in the original lattice setting) and, with more effort and more recently, 
for the Poisson model has been proven rigorously. The references given in Section 5.1 below 
may be used as a starting point into the enormous literature on localization for the Anderson 
model. For the mult i- dimensional Poisson model localization at low energy has been proven, 
separately for the case of positive and negative single-site potentials, in [11] and [12]. 

A localization proof for the RDM provides additional challenges and does not follow from 
the methods alone which have led to proofs for the Anderson and Poisson models. The main 
reason for this is that the RDM does not have any obvious monotonicity properties with 
respect to the random parameters. Such properties are frequently used in an essential way in 
the theory of the Anderson model and, to some extend, can also be exploited for the Poisson 
model. 

This becomes apparent most immediately when attempting to characterize the bottom of 
the spectrum for these models. For the Anderson model with sign-definite potential q this 
corresponds to choosing all couplings minimal (if q is positive) or maximal (if q is negative), 
respectively. For the Poisson model the bottom of the spectrum is if g is positive (due to 
large regions devoid of any Poisson points) and — oo if g is negative (due to dense clusters of 
Poisson points). 

Identifying a mechanism which characterizes the bottom of the spectrum, a crucial prelim- 
inary step towards the localization question, poses a non-trivial challenge for the RDM. It 
is not at all obvious which configurations u = (ojj) of the displacements should characterize 
minimal energy. 

Far more than a characterization of the spectral minimum is needed for a proof of localiza- 
tion. From the theory of the Anderson model it is well known that sufficient ingredients are 
smallness of the integrated density of states (IDS) at the bottom of the spectrum (e.g. in the 
form of Lifshitz tails), and sensitivity of the spectrum to the random parameters (for example 
in the form of spectral averaging or Wegner estimates). The usual approaches to verifying 
both of these ingredients make heavy use of monotonicity properties as well, providing further 
obstacles to a localization proof for the RDM, where such properties are not apparent. 

There are two previous works in which localization properties of modified versions of the 
RDM (1), (2) have been shown. In [19], Klopp considered a semiclassical version — /i^A -|- V^ 
of (1) and was able to show the existence of a localized region in the spectrum for sufficiently 
small value of the semiclassical parameter h. In this regime neither Lifshitz tails nor an 
exact characterization of the spectral minimum are needed for the localization proof. More 
recently, Ghribi and Klopp [14] considered an RDM of the form (1) with an additional peri- 
odic background potential. For a generic non-zero choice of the latter and sufficiently small 
displacement vectors Ui, they use first order perturbation arguments to recover monotonicity 
properties which lead to Lifshitz tails as well as a Wegner estimate, and thus localization. 

Our goal here is to prove localization for the model (1), (2) without working in the semi- 
classical regime or modifying the background. This means that we can not hope to reveal any 
monotonicity properties by exclusively using perturbative arguments, at least not easily and 
not with first order perturbation theory. 
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We will make central use of symmetry properties. In particular, we will assume throughout 
that the real-valued single-site potential q is reflection symmetric in each coordinate, i.e. 



g(. . . ,Xfc_i, —Xk, 



• ^^fc+l) • • •) —?(•••; ^k-lj ^ki ^k+li ■ ■ -ji 



.d. 



(3) 



Due to the use of some higher order perturbation theory we will need some smoothness of 
q and for convenience assume that 

q e C°^(E.'^), (4) 

which could be substantially weakened. We also assume that 

supp q C [— r, r] for some r < 1/4. (5) 

The latter is best understood in conjunction with the following assumption on the dis- 
placement parameters u = {ui)i^x'''- They are i.i.d. M'^- valued random variables, distributed 
according to a measure /i satisfying 



C C SUpp/i C [-dmax^d^ax]'^, 



(6) 



where d. 



max — 2 ~ ^' Here C := {{±dmax, ■ ■ ■ , ^dmax)} dcuotes the 2^ corners of the cube 
'^. For the proof of our main result. Theorem 1.2 below, we will need some 



max: ('"max 



[-d 

smoothness of the distribution of /i. A convenient assumption is that /i has a density p, i.e. 

(7) 



li{B) = Jg p((^o) duo for all Borel subsets B of [—dmax, dmax]'^^ and that 



P fc t-' (,[ 0,maxi(^max\ )■ 

However, as will be discussed in Section 5, our proofs work for considerably more general 
distributions. In particular, our proof of Theorem 1.2 below only uses that fi has a C^-density 
in a neighborhood of the corners C and can be arbitrary away from this neighborhood. 

By (5) and (6), the centers of the "bubbles" q{x — i — cuj) can move all the way into the 
corners i + C oi i + [—dmax, dmaxY, and the supports of the bumps stay, up to touching, 
mutually disjoint. See Figure 1 for a typical configuration. Note that only r < 1/2 would 
be needed to do this in a non-trivial way (give the bubble space to move), but that we will 
need r < 1/4 for technical reasons to make use of results in [2] and [21], see the discussion 
preceding Corollary 3.4 below. 



Figure 1. A typical configuration 
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The random operator if^ is ergodic with respect to shifts in if- and thus, by the general 
theory of ergodic operators (see e.g. [25]), its spectrum is almost surely deterministic, i.e. there 
is a S C M such that 

(t{Hi^) = S almost surely. 

Under the above assumptions it was shown in [1] that among all configurations u a config- 
uration with lowest spectral minimum is given by u* = {u* 



iJidi 



where 



^; = ((-1)^^^. 






(^i,...,^rf)GZ^ 



(8) 



As discussed in [1], under the condition (6) this also characterizes the minimum of the 
almost sure spectrum of the RDM, i.e. 

Eo:=infS = infa(i7^*)- 

In the configuration u* the single site potentials in V^. form densest possible clusters where 
2°^ neighboring bumps move into adjacent corners of their unit cells, see Figure 2. 




Figure 2. The minimizing configuration 



Thus, for the set of assumptions listed above, we have an answer to the preliminary question 
of characterizing Eq = inf S. For a proof of localization we will need much more information. 
In particular, we will need quantitative bounds on the probability that other configurations 
have spectral minimum close to Eq. This will require one more condition on q. 

To state this condition, we need to introduce Neumann operators where a single bump is 
placed into a unit cell at varying position. For this, and for later, set 

Ar = Ar{0) = y--,-j , Ar{x) = Ar + X, Xx = XAiix), 

the characteristic function of the unit cube centered at x. 

rOr O t [ (J'maxi (^maxj i6t 

Hl{a) = -A + q{x-a) on L^{Ai) (9) 

with Neumann boundary condition. Finally, let 

be the lowest eigenvalue of Hj!^ [a). Note that, by symmetry of q, Eq^o) is symmetric with 
respect to all the coordinate hyperplanes Ofc = 0, A; = 1, . . . , d. 
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In [1] the following alternative was established (only requiring boundedness of g, not smooth- 
ness) : 

(i) Either -£"0(0) is strictly maximized at a = and strictly minimized in the corners C of 
\-(] d ¥ 

[ ^maxi ^max\ ; 

(ii) or Eq^o) is identically zero, and the ground state of Hj^ (a) is constant near the bound- 
ary of Ai. 

From this and the symmetries of q it is easy to see, e.g. [1], that u* given by (8) is a 
spectrally minimizing configuration and that the almost sure spectral minimum Eq of Hi_j is 
given by the minimum value of Eo^a), i.e. its value at a G C. 

While it is possible to construct non-vanishing q where (ii) holds, e.g. our remarks in Sec- 
tion 5.2, alternative (i) is the generic case. It holds if g 7^ is sign-definite (since then the 
ground state energy of the Neumann Laplacian —A^ must be shifted up or down), but also 
for generic sign-indefinite q. 

We are now able to state our main result on localization. Here, for a self-adjoint operator H 
and Borel function g we define g{H) by the functional calculus. In particular, Xi{H) denotes 
the spectral projection onto / for H. 

Theorem 1.2. Assume that d > 2, u and q satisfy (3) to (7), and that Eq^q) does not vanish 
identically in a e [-dmax.dmaxY ■ 

Then there exists 5 > such that H^ almost surely has pure point spectrum in I = [Eq, Eo+6] 
with exponentially decaying eigenf unctions. Moreover, H^^ is dynamically localized in I, in the 
sense that for every C < 1? there exist C < 00 such that 

E ( sup \\Xxg{H^)xi{H^)Xy\\l] < Ce-l^-^l' (10) 

Vl9l<i J 

for all x,y E Tf' . Here, the supremum is taken over all Borel functions (7 : M — ;■ C which 
satisfy \g\<\ pointwise. 

Note here that dynamical localization in physical sense is a special case of (10), choosing 
(yf(if ) = e~**^ and taking the supremum over t G M. 

The proof of Theorem 1.2 proceeds via multiscale analysis and our main task will be to 
establish the two main ingredients into the multiscale analysis, i.e. a smallness bound on the 
probability that finite volume restrictions of -f/^^ have low lying eigenvalues (related to Lifshitz 
tails of the integrated density of states) and a Wegner estimate. 

It is in the proof of these two ingredients where new ideas are needed. In both proofs we will 
use that alternative (i) can be strengthened if some smoothness is assumed for q: in this case 
it can be shown that the first partial derivatives of Eo{a) are non-zero in all directions away 
from the symmetry planes. In particular, at its strict minima in the corners C, the function 
Eo{a) is not fiat as it has non-vanishing gradient. The proof of this result, given in Section 2, 
starts from a second-order perturbation theory formula. 

This can be considered as the crucial monotonicity property which makes the localization 
proof work. In the spectrally minimizing configuration u* of the RDM all bubbles sit in 
corners. Non- vanishing of the gradient of EQ{a) in the corners will allow to gain quantitative 
control on how close the ground state energies for other configurations are to Eq. 

First, in Section 3, this will lead to a Lifshitz tail bound with a proof which is based on 
an argument in [21]. This also uses the fact that under alternative (i) the configuration cu* 
is, up to translation, the unique periodic configuration with spectral minimum Eq, a result 
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established in [2], where it is also shown that this is noi true for d = 1. In [21] it was shown that 
this uniqueness result leads to a Lifshitz tail bound for the IDS near Eq if the distribution of 
the Ui is discrete and contains all corners C. In fact, for technical reasons the bound obtained 
there is weaker than the "classical" Lifshitz tail, but it is strong enough for a localization proof 
(if one also has a Wegner estimate). In Section 3 we will show how non- vanishing of VEo{a) 
in the corners can be used to get same bound on the IDS under the assumptions considered 
here, i.e. with distribution satisfying (6). 

What allows us to push the localization argument through is that under assumption (6) we 
can also prove a Wegner estimate, see Section 4. This will again use that VEo^a) 7^ at the 
corners, which will provide us with a measure for how much the ground state energy of finite 
volume restrictions of H^^ is pushed upwards if the bubbles q{x — i — Ui) move away from the 
corners towards the center i of Ai(z). This is technically implemented in Proposition 4.3 below 
in form of positivity of the derivative of H^^ with respect to a suitable vector field, whose proof is 
close to an argument previously used in [22] for operators with random magnetic fields. Based 
on this result, we are able to prove a Wegner estimate by modifying an argument developed in 
[5] and [15] in the context of the Anderson model, which itself is a modification of the original 
argument due to Wegner. In particular, we obtain the correct (linear) volume dependence 
and can conclude Holder continuity of the IDS as a by-product. 

We finally mention that the connection between monotonicity properties of if^^ (a) and the 
RDM is made through surprisingly simple Neumann bracketing arguments, making crucial use 
of the variational characterization of the lowest eigenvalue as the minimum of the quadratic 
form. Thus, we are able to deduce monotonicity properties of a model with infinitely many 
parameters from a one-parameter model. This trick, employed in Sections 3 and 4, was 
previously used in [20] to study the Anderson model with sign-indefinite single-site potential. 
It is the main reason why symmetry of the single-site potential is important for us. 

Our concluding Section 5 serves two purposes. First, in Section 5.1, we briefly discuss how 
the Lifshitz tail bound and Wegner estimate obtained here lead to a proof of spectral and dy- 
namical localization via multi-scale analysis. In Section 5.2 we mention some generalizations, 
related results and open problems. 

2. Bubbles tend to the corners 

In this section we will prove a property of the derivatives of Eo{a), the ground state of 
the Neumann operators H^ (a) defined in (9), which is crucial for all our later arguments. 
The main result of [1], i.e. that under the generic alternative (i) the function Eo{a) is strictly 
minimized in the corners C, may be dubbed as "bubbles tend to the corners". The seemingly 
small but important improvement to be shown here means that "as bubbles tend to the corners 
the rate of change of Eo{a) does not vanish". 

In fact, in [1] two methods were developed to prove that the minimium of Eq^o) is in C. 
One method, used to prove Theorem 1.3 in [1], applied to Neumann operators on rectangular 
domains as considered here and is mostly based on exploiting symmetries of the domain and 
the potential. A second, very different, method was developed in [1] to show the phenomenon 
that "bubbles tend to the boundary" for Neumann operators on general smooth domains and 
with smooth potentials, see the proof of Theorem 1.4 in [1]. What we do here is to apply 
the second method in the setting of rectangular domains. That "smooth methods" apply to 
rectangles is possible due to our use of Neumann boundary conditions, which allow to get 
smooth extensions of the ground state eigenfunction by reflection. The specific geometry of 
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rectangles allows us to push the method used in Theorem 1.4 of [1] further and to obtain the 
non- vanishing of derivatives of Eo{a) in the corners. 

2.1. Basic smoothness properties. Throughout this section we fix an open rectangular 
parallelepiped D G M."^ and consider operators on L'^{D). Let g be a real- valued smooth 
function with closed support in D, i.e. q G C°°{D) such that q{x) = for a: in a neighborhood 
of dD. Consider the quadratic form 

||Vm||2 + {u,qu) 

where we use the symbol 



{u,v) = / u{x)v{x)dx . 
Jd 

This quadratic form is a closed form on H^{D) and defines a unique self-adjoint operator 

H := -Ayv + g 

where A^ is the Neumann Laplacian. The eigenvalues of this operator are discrete, have finite 
degeneracy and tend to infinity. It will be important for us that the eigenfunctions are regular 
up to and including the boundary, i.e., C°°{D). We state this fact as a lemma. 

Lemma 2.1. The eigenfunctions of the operator H can be extended to a neighborhood of D 
where they are infinitely often differentiable. 

Proof. By reflecting the potential and eigenfunctions repeatedly across the boundary one 
obtains a generalized eigenfunction with the same eigenvalue on the whole space W^. Since 
the potential is C°° it follows, by elliptic regularity, that the generalized eigenfunctions are in 

C°°(M'='). D 

We set for aeW^ 

<la{x) = q{x - a) . 

The set of points a for which the support of qa is a subset of D is denoted by G. This is an 
open rectangle as well. Let 

Ha '■= — Aat + qa, 

which is short for the notation H^{a) used elsewhere in this paper, and denote its eigenvalues 
in increasing order and counted with multiplicity by 

EJa) , n = 0,1,2,... 



/, J., , 



and the corresponding real-valued normalized eigenfunctions as 

Un{x,a) . 

In what follows we often denote the extension of the function to a larger set by the same 
symbol. 

Lemma 2.2 (Differentiability of eigenvalue and eigenfunction). The eigenvalue EQ{a) 
as well as the eigenfunction uo{-,a) are (as an L?'[D) -valued function) infinitely often differ- 
entiable in a neighborhood of G. In particular these two functions are in G°°{G). 
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Proof. Continuity of the eigenvalues was proved in Lemma 2.1 in [1]. Since the ground state 
eigenvalue is not degenerate, Eo{a) is strictly separated from the rest of the spectrum locally 
uniformly in a. Hence, there is a circle C with center EQ{a) which is strictly separated from 
the rest of the spectrum, locally uniformly in a. Using the formula, which holds for all z in 
the intersection of the resolvent sets of Ha and Ha', 

{Ha - Z)-' - {Ha' - Z)-' = {Ha - z)"' [g,, - Qa] {Ha' - z)-' , (11) 

we learn that the projection onto the eigenfunction uo{x,a), 



^{{Ha 

^TT? Jc 



Po{a) = —- (i(H.- zTHz 



is also continuous. Repeated use of (11) shows that the projection Po{a) can be arbitrarily 
often differentiated with respect to a. Note that formula (11) holds for all positions of the 
potential, in particular, the support of the potential does not have to be in D. The eigenvalue 
equation 

{Ha - z)-'Po{a) = {Eo{a) - z)-'Po{a) 
now shows that EQ{a), in turn, can also be differentiated as often as we please. D 

Given the above smoothness properties and the non-degeneracy of Eo{a) we can derive the 
following first and second order perturbation formulas for -E'o(o)) 

da^Eo{a) = -(uo, {d^^qa)uo), (12) 

and 

dlE,{a) = {no, {dl^qa)uo) - 2 J^ ^"°'f ^.'^J"'^' , (13) 

for j = 1, . . . , d. This is done in complete analogy to the better known case of non-degenerate 
eigenvalues of operators of the form A+\B, with (12) corresponding to the Feynman-Hellmann 
formula, see e.g. Section 2.3 of [1]. 

2.2. Second order perturbation theory. Making notations more explicit, we now write 

Z} = (ai,/3i)xD', 

where 

^' = («2,/32)x...x(arf,/3,)cM'^"^ 

Then G = I x G", with an open interval / = (— 5i, ^2), 61 > 0, ^2 > 0, and an open rectangular 
parallelepiped G' C M'^"^ 

We will consider the dependence of £'o(ai, . . . , a^) on the first variable ai G / at fixed value 
of (02, . . . , ad) G G'. Throughout this subsection we will abuse notation and write g^^ := g^, 
Hai := — Atv + Qai, as well as -En(«i) and m„(x, ai) for its eigenvalues and eigenfunctions. 

For simplicity of presentation we focus on the variable ai. Clearly, the arguments below 
also apply to the dependence of Eq on each of the other variables 02, . . . , a^. 

The following lemma provides us with a differential equation for the eigenvalue EQ{ai). 

Lemma 2.3. The ground state energy satisfies the equation 

e;; - 4(.„, d,u,)E;, = -2^, ?^f^ . (u) 

'—^ hii. — rL/n 
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Here B{u, v) := {u, An) — (An, v), E'q and Eq refer to ai-derivatives of Eq, and diUo refers 
to the spatial derivative in xi-direction. 

Proof. The proof is essentially the same as the proof of Lemma 5.1 in [1], Due to the fact that 
we work with the partial derivative in a fixed coordinate direction instead of the gradient and 
that D has flat faces, the argument here is simpler than in [1] by not having to include a term 
related to the curvature of the boundary of D. 

Starting with (13), which is a "partial" version of (25) in [1], one follows the argument there 
and arrives at a partial version of (27) in [1], 

Eq = 2'^B{uk,diUo){uk,diUo) -4:{uo,{diqa^)uo){uo,diUo) 

k 

^ Ek-Eo ' ^ ' 

As in [1] the first term term can be rewritten as 
^5(nfc,9ino)(nfc,9ino) = -(^iUq, (-A + g^, - Eo)9ino) + ^(^fc - ^o)|(^fc,5i%)P • (16) 

k k 

Since diUo is in the form domain, 

J2iEk - Eo)\{uk, diUo)\^ + EoWdimW^ = \\{H - Eoy/^dmof + ^oll^i^oir 

[\VdiUo\^ + QaAdiUoY] dx 



'D 

by Kato's form representation theorem [17]. We find from (16) that 

y2B{uk,diUo){uk,diUQ) = [(^lUo, A^iuo) + llVSiUof] = / V ■ {diUoV diUo) dx , 

k ■'D 

where Green's identity was applied. This equals 

E 



{diUo)djdiUodS — / {diUQ)djdiU(jdS 



Sj JT, 



(17) 



where Sj and Tj are the faces of D perpendicular to the j direction. Since djUQ = on Sj 
and Tj we find that the sum can be restricted to the indices j ^ 1. Since djdiUQ = didjU^ 
and djUQ = on 5*^ and Tj, the expression (17) vanishes. After substituting this and the first 
order perturbation formula (12) for j = 1 into (15) we arrive at (14). 

D 

Lemma 2.4. Assume that the right hand side of (14) vanishes for some ai^ G /. Then 
Eo{aifi) = Eq^giq) = and djdiUQ^-, ai^) = on Sj and Tj for all j = 1, . . . , d. 
Moreover, Mo(-,ai,o) is constant on the right and left faces Si and Ti of D. 

It will follow in Theorem 2.5 below that Mo(", «i.o) takes the same value on 5*1 and Ti. 

Proof. The vanishing of the right side of (14) means that, for k ^ 0, 

(Awfc, diuo) - {uk, AdiUo) = . (18) 
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Using the eigenvalue equation —Au^ + Qai o^fc = -^fc^fc "we can transform this into 

(Efc - Eo){uk,diUo) = -(m/c, idiqa^o)uo) ■ (19) 

Now pick any function / G C°°{D) and use the previous identity to write 

^{Ek - Eo){uk,diUo){f,Uk) = -^(wfc, (<9igai,o)Mo)(/,Mfc) + {uo,diqa,,oUo)if,Uo) . (20) 

k k 

Since the functions Uk-, A; = 0, 1, 2, . . . form an orthonormal basis for L'^{D) and /, duo are in 
the form domain we use once more Kato's representation theorem to obtain 

'^{Ek-Eo){uk,diUo){f,Uk) = (V/, V((9iMo)) + {f,qai,Auo) - Eo{f,diUo) . 

k 

Plugging the latter two identities back into (20) we arrive at 

(V/, V(9iMo)) + (/,5l(gai,oMo)) -Eo{f,diUo) = (Mo, (5igai,o)Mo)(/,Mo) • (21) 

Now pick / G C^{D) so that we can integrate by parts without boundary terms and get 

(/,(9i(-Amo + gai,o^*0 - EqUo)) = (mo, {d\qa^,o)uQ) {f ,Uq) . (22) 

Since the left side vanishes for all / G C^{D) and the latter are dense in L'^{D), we must 
necessarily have that 

{uq, {diQa^Juo) = , 

i.e., ^o(«i,o) = by (12). Thus, we have for all / G C°°(L>), 

(V/, V(9ino)) + (/,9i(g,,,„no)) - i?o(/,5i%) = . (23) 

Doing the same integration by parts as above, this time with boundary terms, yields 



E 



J 



{d,d,uo)fdS- / id,d,uo)fdS 

Si JT-i 



(24) 



for all / G C°^{D). This means that 

djdiUQ = (25) 

pointwise on Sj and Tj for all j = 1, . . . , d. In particular, dfuo{x, ao,i) = on 5"! and Ti. Since 
the potential vanishes on the faces 5*1 and Ti we find that Uq satisfies the equation 

-A'uo = -J2Y^ = Eo{ao)uo (26) 

on 5*1 and Ti. The function uq is smooth up to and including the boundary of -D, in particular 
it is a smooth function on the faces of D. Consider the d — 2 dimensional 'edge' where 5*1 and 
5*2, say, meet. The gradient of uq at this intersection must be of the form 

Vuo = {0,0,dsUo,---ddUo) . (27) 

Hence, Uq is an eigenfunction of —A' on 5*1 and Ti with a Neumann condition on the boundary. 
Since uq has a fixed sign, i?o(a'i,o) must be the lowest eigenvalue of the Neumann Laplacian 
on 5*1 and Ti and hence 

^o(«i,o) = , (28) 

and Mo is constant on 5*1 as well as on Ti. D 
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Theorem 2.5. Assume that the right side of (I4) vanishes for some a^'^' E G. Then Eo{a) = 
identically in G and for every a E G the ground state Uo{x,a) of Ha is constant near the 
boundary of D. 

Proof. Write a^°^ = (ai,o, 02,0, • • • , a^.o) and consider g^^, -E'o(ai) and Uo{x,ai) as functions 
of the first coordinate of a = (ai, 02, . . . , a^) only, setting (02, . . . , a^) = (02,0, • • • , 0^,0) in 
Lemmas 2.3 and 2.4 above. 

We have ai^o ^ (^<^i;<^2) and thus the faces Si and Ti of D are in the complement of the 
support of gai,o- 

Pick a point on 5*1. We may call this point the origin. By the reflection argument already 
used in the proof of Lemma 2.1 the function uo{x,aifi) is harmonic in a neighborhood A^ of 
0. Thus, it has a convergent power series expansion 

:c,a;" (29) 



E' 



which can also be written as 



E' 



C!0 



>^2^C(fc,«,)X- =}_^xtP,{x'), (30) 

fc=0 a' k=0 

where x' = {x2, ■ ■ ■ , Xd)- Since Uq is harmonic we find the recursion 

A'Pfe(x') + {k + l){k + 2)Pk+2{x') = , A; = 0, 1, 2, . . . . (31) 

By Lemma 2.4 we know that no(0,x', Oi^o) is constant and we find that Po{x') is constant, 
too, and thus Pk{x') = for fc > 2, even. Since 9i9jMo(0,x') = for j = 2, ... ,d we find 
that djPi{x') = for j = 2, . . . , d and hence Pi(x') is constant and thus on account of (31), 
Pk{x') = for all k > 2. Thus, uo{x, ao,i) = c + dxi near 0, where c, d are constants. However, 
as Mo(x, ao,i) satisfies Neumann boundary conditions, we must have d = 0. This implies that 
uo{x, ao,i) = c on (ai, ai + 61 + ai,o) x D', where ga^ ^ = 0. 

In the same way, now starting with a point on the face Ti it is shown that uo{x,aifi) is 
constant on (/3i — ^2 + ai,o, A) x D'. 

If we translate the potential ^^1,0(3;) — ^ Qai^oi^ ~ l*^! ~ '^1,0)) = Qi^ ~ 0,1), the translated 
function v{x) = uo{x — (ai — aifl), aifl) is an eigenfunction of the Schrodinger equation Ha^v = 
Eo{aifi)v and since v has fixed positive sign, it is the ground state. Thus, the eigenvalue 
Eo{ai) = identically in ai G /. 

Note that, by assumption in Theorem 2.5, the above argument can be applied to Eo{ai, . . . , ad) 
as a function of each one of its variables with the other coordinates fixed. Thus, it vanishes 
identically in each variable and therefore identically on G. Also, it is found that UQ{x,a^'^') 
is constant in a neighborhood of each face Sj, Tj, j = 1, . . . ,d. As the union of these neigh- 
borhoods is connected, it follows by unique continuation that Uo{x,a^^^) is constant near the 
boundary of G. This is then also true for the ground state uo{x, a) = uo{x — (a — oq), ao) with 
arbitrary a E G. D 

For what follows, it will be convenient to assume that (ai,/3i) = (— s, s) for some s > 0, 
supp q C [— r, r] x D' for some < r < s, and that q is reflection-symmetric with respect to 
the first variable, i.e. q{xi, x') = g(— Xi, x') for all Xi E (— s, s) and x' E D' . As a consequence, 
the function Oi — )■ -Eo('2i) is symmetric about oi = 0, in particular E'q{0) = 0. 

Corollary 2.6. Let q be smooth and be reflecton symmetric with respect to the first variable. 
Assume that -E'o(a) does not vanish identically on G. Then, for each choice 0/(02, • • • , CLd) E G' , 
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the function Eo^ai) = Eo{ai, a2, . . . , a^) satisfies Eq^gi) < for all ai G (0, s — r] and, by 
symmetry, Eq{0) = and E'^lai) > for ai G [— (s — r), 0). 

Proof. For fixed (02, . . . , a^) G G', denote the right hand side of (14) by C(ai). By assumption 
and Theorem 2.5 we know that C(ai) is strictly negative for ai G (— (s — r), s — r). Using the 
integrating factor e^^"^-*, where F{ai) is an antiderivative of — 4(no, 9ino)(ai), equation (14) 
can be written as 

-^ (e^(-)i?^(ai)) = C(aOe^('^^) . 

Integrating this equation using the fact that E'q{0) = yields the result. Note that this 
holds up to the boundary of [—{s — r),s — r] by smoothness of Eq^q) up to the boundary 
(Lemma 2.2). D 

We are ready to state the strengthened version of Theorem 1.3 of [1], which will be central 
to all our later considerations in the localization proof. 

The domain D in Mf^ may be any open rectangular parallelepiped and q a C°°-smooth 
function with closed support contained in D and having the same symmetry hyper-planes 
as D. We define G and Eo{a) for a G G as in Section 2.1. Then G is an open rectangular 
parallelepiped and we may assume that its center is in the origin and G = {—Mi, Mi) x . . . x 
{—Md,Md). Due to the symmetry assumption, Corollary 2.6 applies to each coordinate of 
Eo{a) and yields 

Corollary 2.7. Assume that Eq^q) is not identically zero on G. For all a = (ai, . . . ,ad) G G 
and all i = 1, . . . ,d we have 

( <0, ifai> 0, 
diEo{a) 1=0, if ai = 0, 

[ > 0, if Qi < 0. 

Open Problem: It is an open problem to show that the phenomenon "bubbles tend to the 
corners", i.e. results such as Theorem 1.3 in [1] or Corollary 2.7 above, also appears if the 
rectangular domain D is replaced by more general polyhedra, for example regular n-gons in 
M^. Again this should need a suitable symmetry assumption on the potential q, e.g. spherical 
symmetry. Particularly interesting would be the case of an equilateral triangle, as all other 
results in this paper would be applicable to this case as well. It was shown in [1] for general 
smooth convex domains that the minimizing position lies at the boundary. 

3. LiFSHITZ TAILS 

In proving smallness of the IDS of H^^ near Eq, the almost sure spectral minimum of H^, 
we will rely on a prior result of this kind for the RDM (1), (2) established in [21]. Their result 
requires widely the same assumptions as made above, i.e. (3), (5), (6) and that Eo^a) does 
not vanish identically, but needs in addition that the support of the distribution /i is a finite 
discrete set. This violates our assumption (7) (and every other assumption on /i which works 
for our proof of a Wegner estimate in Section 4). 

However, based on the results from Section 2, one might expect that configurations in which 
all bubbles sit in corners, as dealt with in [21], constitute the worst case scenario for smallness 
of the IDS near Eq. We will make this rigorous here by showing that the IDS of H^^ can be 
bounded from above (up to a multiplicative constant) by the IDS of a modified RDM where 
all bubbles have been moved to the nearest corner. 
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More precisely, for an a G G = [—dmax, dmaxY, let c{a) denote the element of C closest to a. 
If there is more than one such point, any of them can be chosen; for the sake of definiteness, we 
may order the points in C lexicographically and chose the first in the list. For a displacement 

configuration u = {uji)i^id G G we will write c{uj) for the closest corner configuration given 
by {c{uj))i = c{uji), i G Z*^. 

For a non-negative integer L, let A2L+1 = {—L — 1/2, L + 1/2)"^ and H^^l the restriction of 
H;^ to Al with Neumann boundary conditions. Also, let A2^_^^ = Z'^ fl A2L+1. We will prove 

Proposition 3.1. There exists a constant C G (0, 00) such that, in the sense of quadratic 
forms, 

Huj,L — Eq> -^(-f/^c(a;),L — -^o) (32) 

^d 

for all uj G G and all L > 0. 

We start by showing this for L = 0, i.e. on the level of the single-site operators H^ (a). This 
is achieved in the following two lemmas, which separately treat the cases where the bubble is 
close to a corner or not close. 

Recall from Section 1 that EQ{a) denotes the ground state energy of Hj^ (a) for a E G and 
that the almost sure spectral minimum Eo of H^ is given by Eo{a) when a is one of the corners 
C of G. As a consequence of Corollary 2.7, EQ{a) grows linearly in the distance of a from the 
nearest corner, i.e. there exists a constant G G (0, 00) such that, 

^o(a)-^o>^^(a) (33) 

where D{a) = miucgc |a. — c|. 
We will use this to prove 

Lemma 3.2. There exists G > and 6 > such that, if D{a) < 5, then 

HU<^) -Eo>^ «(c) -Eo+\a-c\)). (34) 

Proof. Let a E G and pick c E C such that D(a) = \a — c\. As q is G^, write 

Hlia) -Eo = Hlic) -Eo + q{- ~ a) - q{- - c) 

= Hlic) -Eo + {c-a)- Vg(- - c) + o{\a - c\). 

Using (33), one obtains that 



(35) 



Hlic) -Eo + {c-a)- Vg(- - c) > -\a - c\ + o{\a - c\). 

Hence, there exists p E (0, 1) small such that, for a E E>'^~^ with a = c + pa E G, one has 

Hlic) -Eo-pa- Vg(- - c) > p/2G. 

We recall Lemma 2.1 of [20]: If we suppose A > and A + B > cq > 0, then we have 
A + tB> min(l/2, cq) -{A + t) for t E [0, 1/2], since 

A + tB = {l~t)A + t{A + B)>lA + tco > min(i, co){A + t). 
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Applying this with A = Hj^ {c) — Eq and B = —pa ■ Vg(- — c), we learn that there exists 
Cp = niax(2, 2C/p) > such that for t G [0, 1/2] and for a G S'^^^ with c + pa E G, one has 

Hlic) -Eo- tpa ■ Vg(- - c) > -^ (hUc) - Eo + t) . 

Hence, in view of (35), for |a — c| < p/2 and t = \a — c\/p, one has 

Hlia) -E,>^ (Hlic) -Eo + \a- c\/p) + o(|a - c|). 

Finally, this implies that there exists 6 > such that for \a — c\ < 6, one has (34) and completes 
the proof of Lemma 3.2. D 

Bubbles which are not close to a corner are easier to handle and considered in 

Lemma 3.3. Fix 6 G (0, 1). There exists C = C^ E (0, oo) such that, for D[a) > 6 and all 
c E C, one has 

HZia) -Eo>^ {Hlic) -Eo + \a-c\). 

Proof. Fix 6 G (0, 1). By (33) there exists r/ > such that, for a with D{a) > 6, H^^{a)-EQ > 
rj. Hence, as \q{x — a) — q{x — c)\ > — 2||g||oo, there exists C > 1 such that, for a satisfying 
D{a) > 5 and all c E C, one has 

(C + l){Hl{a) - Eo) - (Hlic) - Eo) > Cr] - 2||g|U >V>^\a-c\. 
Hence, 

H^a) -Eo> -^ (^Hlic) -Eo + ^\a-c\ 

>C(^^iKic)-Eo + \a-c\). 

This proves Lemma 3.3. D 

To complete the proof of Proposition 3.1 we need to extend the result of the previous 
lemmas to general boxes A2l+i- This is done by an argument previously used in the proof of 
Theorem 2.1 in [20] which makes crucial use of properties of Neumann boundary conditions. 

For ip G H^{h.2L+i)i the form domain of H^^l., one has that the restriction of ^ to Ai(z) is 
in H^{Ki{i)) for each i G Ag^^^^ and 

{{H^,L-Eo)^,ij)= Y^ ((-A-Eo + g(--i-Wi)^,^)Ai(i) 

where (■, ■)a_ denotes the standard scalar product in L'^{A). This may also be applied to the 
modified displacement model Hc(uj), 

((^cH,L -^o)^,^) = Yl ((-^-^o + g(- -«-c(a;i))^,^)Ai(i)- 
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Hence, using Lemmas 3.2 and 3.3 on each term in the sum, we obtain that 

where the positive term X]i(l'-^« ~ c{ui)\ip,ip)\-^(^i) was omitted. This completes the proof of 
Proposition 3.1. 

The random displacement model -ffc(aj) = —^ + J2iez^Qi' '~ c(wi)) has i.i.d. displacement 
vectors (c(ci;))jgzd, whose distribution is discrete with support given by the corners C of G. 
Thus, it satisfies the assumptions of Theorem 4.1 in [21]; hence, by the proof of Theorem 1.2 
in [21], in particular (3.2) in the same work, there exist Ci > 0, C2 < oo and /i > 1 such that, 
for all L, 

¥{Hc(uj),L has an eigenvalue less than Eq + Ci/ L^) < C2L'^/i~'^. 

Note that this requires that d > 2. The argument leading to Theorem 4.1 of [21] uses 
crucially the uniqueness (up to translations) of the minimizing configuration of the potentials 
proved in [2], which holds only for d > 2. See Section 5.2 below for a comment on the 
differences for d = 1. The assumption r < 1/4 in (5) was used in [2] for a more technical 
reason (rather than just r < 1/2) and thus also enters our argument here. 

Using Proposition 3.1, we immediately obtain the following finite volume bound on the prob- 
ability for finding low lying eigenvalues. It is this result which enters the proof of localization 
via multiscale analysis. 

Corollary 3.4. There exist Ci > 0, C2 < 00 and fi > 1 such that, for all L, 

'^{Huj,L has an eigenvalue less than Eq + Ci/ L?) < C2L fjT . (36) 

The results of [21] also show that the integrated density of states of i/c^, say N, satisfies 
a Lifshitz tails estimate of the form 

r log I log iV(i?) I ^ 1 

limsup— — -, — < — . 37 

E>Eo 

By Proposition 3.1 we see that A^, the integrated density of states of H^^, satisfies, for E > Eq, 

N{E) <N{Eo + C{E-Eo)). 

Hence, we have proven the following result which is not required in the localization proof but 
stated here for its independent interest. 

Theorem 3.5. The IDS N of Hi^ has a Lifshitz tail of the form 

\og\hgNiE)\ 1 

limsup — — 7— — ^ < — . 

E^Eo log{E-Eo) - 2 

E>Eo 

We expect that the Lifshitz exponent 1/2 is not optimal and should instead be d/2, the 
standard value known from the Anderson model. We can think of two ways in which one could 
try to get this improvement, none of which we know how to make rigorous. One approach 
would be to show that (37) holds with exponent d/2, which would immediately give the same 
in Theorem 3.5. That this should hold is discussed in [21], where the exponent 1/2 is found 
due to one part of the proof which uses an essentially one-dimensional argument. 
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Another way to argue would be to make use of the term J2i \^i ~ c(a;j)|xAi(i) which was 
dropped in the proof of Proposition 3.1. Under our assumptions (6) and (7) this means 
that one would have to show standard Lifshitz tails for an Anderson-type model where the 
unperturbed operator is the random operator Hc(uj)- However, the known methods do not 
work for the irregular background potential appearing here. 

4. A WEGNER ESTIMATE 

Throughout this section, we write Hi{a) = Hj^^,^s{a) for simplicity. Our goal in this section 
if to prove the following Wegner estimate for energies near Eq = mi^^QEQ{a), where EQ{a) = 
ini a{Hi{a)): 

Theorem 4.1. There exists 6 > such that, for any a G (0, 1), there exists Ca > such that, 
for every interval I C [Eq, Eq + 6] and L eN, 

E{tTXi{H^,L))<C^\irL''. (38) 

By Chebychev's inequality this implies that for every interval [E — rjjE + rj] C [Eq, Eq + 6], 

P(dist((7(iJ^,i, E)) < r/) < C;,r/°L^ (39) 

the more classical form of the Wegner estimate used in applications. As a consequence of the 
existence of the integrated density of states (see e.g. [25]) and Theorem 4.1, we also get 

Corollary 4.2. There exists 6 > such that, for any a G (0,1), the integrated density of 
states of H^ is a-Holder continuous in [Eq,Eq + 5\. 

The rest of this section will be devoted to the proof of Theorem 4.1. 
For a function / on G we set 

(5c/)(a):=f|4^-V/(a), 
\c{a) — a\ 

with c{a) denoting the corner closest to a as in Section 3. Thus, dc denotes the directional 
derivative in the direction of the closest corner, where points a with multiple closest corners 
will not play a role in the arguments below (starting from (41) below we introduce a cut-off 
which restricts the values of a relevant for the proof to small neighborhoods of the corners). 
By Corollary 2.7 there exist (5o > and ro > such that 

dcEoia) < -6q for all a e Ar^ := {a e G : \c{a) - a\ < Tq}, (40) 

a neighborhood of C. 

Let rj G C°°(M) such that < ?7 < 1, //(r) = 1 for r < Tq and ?7(r) = for r > 2rQ. Using 
this function as a cut-off, we localize the vector fields associated with dc onto a neighborhood 
of the corners, defining 

{d'J){a):=r^{\c{a)-a\){dJ){a). (41) 

For each 2 G Z*^, we write 

^U^. = d'c,^M- - ^ - ^.) = -^i\c{oo.) - oo.D^^p^^ ■ (Vg)(- - z - c.,). (42) 

li ip E H^{A2L+i), the form domain of H^^, then ipi := '?/'|Ai(i) ^ -^^(^ilO)' ^^e form 
domain of Hi{uji), and, with the usual abuse of notation for the quadratic form, 

{^,H^,Lij)= Yl {i^'-^^H,{ujdi^,), (43) 
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as well as 

Y, (^,9;,,^.,L^)= Yl {i^ud',,.M^^)^^)- (44) 

Proposition 4.3. There exist 5i > and ^2 > such that 

- E (^^K^.H^,L)i^)>64n' (45) 

/or all LeN, andi^e H\A2l+i) with {^j, {H^^l - Eo)i>) < ^sllV'll'- 
In the proof of Theorem 4.1 at the end of this section, Proposition 4.3 provides the crucial 
technical result which will allow us to link our argument with previously known strategies for 
proving Wegner estimates. Thus, we pause here to motivate the origin of Proposition 4.3, 
comparing with the situation of an Anderson-type model 

i7^ = -A + J]M(a;-z) 

with random coupling constants Aj and suitable non-negative single-site potential q. If H^j^ is 
the restriction of H;^ to A2L+1 with appropriate boundary condition, / C M a compact interval, 
and {(pj} all normalized eigenfunctions to eigenvalues of H;^^^ in J, then Proposition 4.3 should 
be considered an analogue to the fact that 

J2 (0,,g(--O0.)>Co>O 

uniformly in j and L, which is a key step in the proof of a Wegner estimate for H^ in J, e.g. 
[5]. 

For this, first note that, in a sense as in (42), dH^j^^/dXi = q{- —i). In the RDM the random 
parameters Ui are vector-valued, which gives us the choice to differentiate in (45) with respect 
to a suitably chosen vector field, for which we take for every Ui the partial derivative in the 
direction of the closest corner. Finally, as indicated by the prime in (45), we need to sum only 
over those i where Ui is close to a corner. This is plausible by the fact that we establish (45) 
only for ip with energy close to Eo, which forces them to have most of their mass in cubes 
where Ui is close to a corner, see Lemma 4.5 below. 

The proof of Proposition 4.3 will be prepared by a series of lemmas. 

Let Pi{a) be the eigenprojection onto the groundstate of Hi{a), Pi{a) := I — Pi{a), and 

E, := ml{a{H,{a)) \ {Eo{a)}) > Eq. 

a&G 

Lemma 4.4. If tp e H^k^L+i) with {^j, {H^^l - ^0)^) < '^211^^11', then 
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Proof. We have 



»GA'2i,+i 






|2 
^t ll ) 



»eA'2L+i 



which yields (46) by the assumption. D 

By the results of [1] we have 

Eo,ro ■= i?f Eo{a) > Eq. 

Define Kl^+i ■= {^ ^ A'a^+i : UJi ^ A^J. 
Lemma 4.5. Ifip G H\A2l+i) with {iP, {H^^l - ^o)^) < h\\M? , then 

E \\n'<jr^^\m'- (47) 

Proof. This follows from the assumption and 



iGA^'i+i 



n 



Lemma 4.6. There exist Ci < oo and C2 < oo such that for (p G L'^{Ai{i)), 

- (y., (9;,i/.(a))(p) > -id',^,Eoia))\\P,iaM' - C,M \\P,cp\\ - C^m^f. (48) 

Proof. Omitting the variable a we have, 

d'^p. = -d[T. = -d',(P-) = P0,P.) + {d',P.)P. 

and 

d'.Hi = d',{E^P, + PiH,7,) 

= {d',Eo)P^ + Eo{d',Pi) + {d'J^i)H,P, + P.id'MPi + P,H0J^{) 
= {^',Eo)P^ + Eo{^',Pi)P, + EoP0^P,) 

+ P0,H,)P, - {d',P,)H,P, - P^R^Pi). 
Thus, forv9GL2(Ai(z)), 

-{if,{d',H,)if) = -{d',Eo){ip,P^f) -2EoRe(P,if,{d',P,)if) 

- (P,^, (d'MP,^) + 2Re{H,{d',P,)v,P,v). 

This implies (48) if we can show that the operators d'^^Hi{a), {d'^^Pi){a) and Hi{a){d'^^Pi){a) 
are bounded in L^(Ai(i)), uniformly in a G G. This is clear for 

^'c,aH^{a) = -^{\c{a) - a|) j^| ~ ^ ■ (Vg)(- -t-a). 
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For Pi{a) we write, as in Section 2, 

Pi(a) = ^^ ([ {Hi{a) - z)-Uz, 

where C circles around Eq^o), and can be chosen locally independent of a and with distance 
to Eq^o) which is bounded below uniformly in a G G (as the distance of EQ{a) and Ei{a) is 
uniformly bounded below). From this we can conclude that 

{K,ama) = ^ £{HM) - zr\&,^,H,{a)){H,{a) - z)"' dz 

is uniformly bounded in a. In this expression one can absorb an additional factor Hi[a) to 
show uniform boundedness of Hi{a){d'^^Pi){a). 

n 

We are now ready to complete the proof of Proposition 4.3: 
Proof. By Lemmas 4.4 and 4.5 we have 



> u\\^ 



> 1 



A" i 

62 S2 



i^^2L+l »e^'2L + l 



El — Eq Eo^ro ~ -^0 

Let C3 := supagfj Wca^oid)]. With Lemma 4.6 we find 

- E {i^AdUHu.,Lm = - E {i^^A^UHM)m 



^^Kl+i 



-Ci E \m\\\Pii^^)i'^\\ - C2 E \\P^i'^^)i'^f 

*^^2L+1 *S^2L+1 

\ 1/2 

Ci„..„/ ^ ||P^(a;,)^,f j -C2 E ll^^(^^)V^^ll'- (49) 



2 



7PA' / jpA' 



Now we use that —d'^^.EQ^Ui) is bounded from below by — C3 if z G f^2L+i ^^^ by Sq if 
z ^ ^2L+i (^be latter means cuj G Aq and thus d'^^EQ^Ui) = dc^^^EQ^oJi), so we can use (40)). 
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Also using the bounds from Lemma 4.4 and 4.5 again, we arrive at 



'6^2L+1 



-C3 



a 



Eo,ro — Eq 2 \ El — Eq 

62 



2" 



' El — Eq 

Choose 62 > such that 



C362 . Ci I 62 , C2S2 So 



-E'o,ro — Eq 2 \ El — Eq Ei — Eq 4 

and 

^2 , S2 5o 



El — Eq Eo^ro — Eq 2 

Then 

This proves (45) with 6i = Sq/A. 



D 



We now prove Theorem 4.1. We follow the approach developed in [15] based on L^ estimates 
of the spectral shift function (see also [5]). The method can be adapted to our model thanks 
to Proposition 4.3. 

For 62 from Proposition 4.3 choose 6 = 62/ 2 and let I C [Eq, Eq + S] be an interval of the 
form I = [E - e,E + e]. 

Let X ^ C°°(]R) be a real- valued function such that x(x) = — 1 for x < —e; xi^) = for 

X > e; x' ^ ^1 and ||x'l|oo < l/^^- 

By our assumption on I and the Gohberg-Krein formula, see e.g. Proposition 2 in [23], 
Proposition 4.3 implies that 



tr 



E d',^Jx{H^,L -E + t)] \=ti ix'{H^,L-E + t)l- Y. {&,^^M. 



>5iti {x\H^,L-E + t)). 
Then, as supp x' C [—e, e] and as x' ^ 0, one has 

r2e 

EitT xiiH^,L)) < E(tr / x'{H^,L-E + t)dt) 

^ tH f^E{tr {-d',JxiH^,L-E + t)]))dt. (50) 

In the above expectation we want to write the integration with respect to Ui over G as 
a sum of integrals over the intersection of G with each one of the 2*^ orthants, using polar 
coordinates with respect to the corners c G C in each orthant. For this we represent a & G 
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by (r, 6', c{a)) G (0, oo) x E>'^~^ x C, where c{a) again denotes the corner closest to a and (r, 6) 
polar coordinates of a — c{a). For a function / supported near the corners C this means that 

f{a)p{a)da = \2 / f {rO - c) p{re - cy^'^ dr dO , 

Jg 7^ Af-^ io 

with S^~^ denoting the intersection of S'^"^ with the orthant containing — c. 
With a = a;,- this leads to 



E(tr (-9;,Jx(^.,L-i5 + t)])) 

= E,/tr J2 [^ [ " dr[x{H^,L-E + t)]r]{r)p{r9-cy-'drd9\ (51) 



where Ej denotes the expectation with respect to the random variables {u}j)j^i. Here a sign- 
change is due to the fact that dr acts in the direction opposite to dc^ui,- By integration by 
parts, we have 

r2ro 

/ dr [x{H^,L -E + t)] r]{r)p{re - c^-^dr 
Jo 

r2ro 

= - / [x{H^,L -E + t)- x{H^i^\L -E + t)] dMr)p{re - c)r''-^)dr, (52) 
Jo 

where u^^^ is the random variable such that a; = tUj for j ^ i, and w^ = c(uJi). Note 
that the second term in (52) actually integrates out to zero, as r]{r)r'^~^ vanishes at both 
endpoints. But we include this term in the integral to be able to make use of bounds on 
X{H^,L -E + t)- x{H^ir).L -E + t). 
Now (52) implies that 



2ro 

tr ( / dr\x{H^^L-E + t)]T]{r)p{re-c)r'^-^dr 

-2rQ 



^~r° (_^e(- ;i^.,L,i^.(.),L)x'(A)rfA') dr[v{r)p{re-c)r''-']dr (53) 

This uses the spectral shift function ^(A ; H^^, H^(i) j^) for the pair {H^^, H^n) j^) which is 
defined so that 

/ <^'(A)e(A ; iJ^,L, i^.w,L)rfA = tr {'^{H^,l) - ¥^{H^i^),L))- 

Jm. 

for all if G C^(M). The invariance principle for the spectral shift function (see e.g. [26] or [3]) 
states that the spectral shift function ,^(A ; H^^l, H^{i)^i^) can be written as 

^(A ]H^^L^H^(i)j} = -^{gk{X) ]gk{H^^L),gk{H^ii)j)). 

Here we define 5'fc(A) = (A + M)~^ and M is picked such that w-a.s., inf^,. K;(x) > —M + 1. 

By definition H^^^l — H^ii)^^^ = q{- ~i ~ ^i) — q{- — i — c{ui)) which is a bounded, compactly 
supported potential. In Section 5 of [15], it is proved that (actually for more general operators), 
a k > pd/2 + 1 and p > I, the operator gk^H^^L) — gk{H^(i)j^) is super trace class of order 1/p, 
i.e. its singular values to the power 1/p are summable; the p-th power of this sum is denoted 
by II ■ ||i/p. Moreover, \\gk{H^,L) - gk{H^(,)^^)\\i/p < Cq for Cq > independent of u and L. 
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Using a simple change of variables and the bound ||^(-; A, -B)||lp < ||A — i?||-^(p proven in [5] 
we find 



Jr 

< CCo. (54) 

As llx'lloo < 1/^ and suppx' C [— £^,e], (54) and the Holder inequality imply that, for any 
q e (1, +CX)), there exists Cg > such that 



sup 



/"e(-;^-,L,^.«,L)x'(A)t/A 

Jr 



<a,e'/^-\ 



As dr[ri{r)p{r9 — c)r'^ ^] is a bounded, compactly supported function uniformly in 9 and c, 
(53) then implies that 

-2ro 



sup 



tr ( / dr[xiH^,L -E + t)]7]{r)p{re - cy~^dr \ 



<a,e^/''^\ 



Plugging this into (51) and then into (50), we get that, for any q G (1,+cx)), there exists 



Cq > such that 

E(tr Xi{H^,l)) < Cq \Kl+i\^'^'-' [ ' dt < Cqs'/'^L''. 

J~2e 

This completes the proof of Theorem 4.1. D 

5. Concluding remarks 

5.1. Remarks on multi-scale analysis. It is well known to experts in the field that Lifshitz 
tails of the IDS, more precisely, a bound as in Corollary 3.4, combined with a Wegner estimate 
such as Theorem 4.1 lead to a proof of spectral and dynamical localization via multi-scale 
analysis (MSA). For the sake of reaching a broader audience we include some discussion of 
the strategy and additional tools which are behind this. 

Two very convenient references for this discussion are [25] and [18], which both make a 
point of thoroughly discussing detailed input assumptions which make MSA work and thus 
apply to a wide range of models, including ours. We follow the survey article [18] here, as the 
results presented there are based on [13], where it was shown that MSA leads to dynamical 
localization in the form (10), the strongest result which has been obtained via MSA. 

The required assumptions singled out in [18] can be divided into deterministic and proba- 
bilistic assumptions. The deterministic assumptions listed in the following hold for classes of 
Schrodinger operators much larger than what we require here. For discussion and references 
on their proofs see [18] or [25]. We use the same acronyms as [18]. 

(i) A property leading to the existence of suitable generalized eigenfunction expansions 
(SGEE). 

(ii) A Simon-Lieb inequaltiy (SLl) relating resolvents at different scales. 

(iii) An eigenfunction decay inequality (EDI) providing estimates for generalized eigenfunc- 
tions in terms of finite volume resolvents. 
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The required probabilistic properties are 

(iv) Z'^-ergodicity, which could be relaxed as discussed in [18], but is clearly satisfied for 
our model. 

(v) Independence at distance (IAD), trivially satisfied in our model due to the non-overlap 
assumption for the single-site terms. 

(vi) A Wegner estimate of the form (38), (39). In [18] the stronger form E(tr XiiH^^ l)) < 
Co,|/|L'^ of the Wegner estimate (linear in the interval length) is used as an assumption, but 
all arguments can be modified to only require the slightly weaker (38), see Remark 4.6 in [18]. 

For random Schrodinger operators satisfying all these assumptions it was shown in [13] that 
a certain "suitability" -property of the finite volume resolvent, sometimes referred to as an 
initial length estimate, implies localization in the form claimed in Theorem 1.2. To establish 
this initial length estimate, the Lifshitz-tail bound from Corollay 3.4 is used. The argument 
behind this is well known, with details found, for example, in [25]. Thus, we only outline the 
two main steps: 

First, by a Neumann bracketing argument one deduces the following from (36): For any 
^ > and (3 G (0, 1) there exists L2 = -^2(/3, such that 

F{H^^L has an eigenvalue less than Eq + L^~^) < L~^ 

for all L > L2. The main difference to (36) is that one trades in a larger distance of eigenvalues 
to Eq for less, but still sufficient, probability. 

Second, for E ^ I := [Eq,Eq + L'^"^/2] one may now use a Combes-Thomas estimate 
(which holds for very general semi-bounded Schrodinger operators, and thus certainly in our 
setting) to turn this into an initial length estimate such as the suitability property required 
in Theorem 5.4 of [18]. This allows to start the MSA machine which leads to all the results 
stated in Theorem 1.2. 

The quadruple MSA needed to prove Theorem 5.4 in [18] is carried out in [13]. An equally 
self-contained but somewhat less refined MSA scheme is provided in [25]. Here MSA is iterated 
twice, obtaining pure point spectrum with exponentially decaying eigenfunctions in the first 
run-through and using a second MSA argument to prove dynamical localization in the form 

E I sup \\\X\PgiH^)xiiH^)Xv\\] < 00 (55) 

Vl.l<i J 

for all p > on an interval / = [Eq, Eq + S{p)] with 6{p) > depending on p. This is much 
weaker than (10), but we still consider the presentation in [25] as a very accessible introduction 
into the mathematics of MSA for non-experts. 

5.2. Related results and problems, (i) The Wegner estimate and the Lifshitz tails, and 
therefore our main result Theorem 1.2, hold under weaker assumptions on the distribution 
fi of the displacements. E.g. the proof as written in Section 4 only requires that /i has a 
C^-density p near the corners C. This is made possible through the introduction of the cut-off 
rj supported near the corners in (41). 

It is evident from (53) that the C^ condition for p is only needed in the radial direction 
with respect to the corners. In fact, similar to [14] we could allow distributions supported on 
a suitable submanifold. Examples would be the uniform distribution supported on a "cross" 
in d = 2 or, in general dimension, uniform distribution supported on the boundary dG of 

'-^ I, (^maxy (^max ) ■ 
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(ii) Under our assumptions and in d >2, among all periodic configurations cu E G , w* as 
defined in (8) is, up to translations, the unique minimizer in the sense that inf a{Hi^*) = inf S. 
This was proven in [2] and enters the argument in [21] leading to Corollary 3.4 and the Lifshitz 
tail estimate Theorem 3.5 above. 

As also shown in [2], in d = 1 there are many periodic minimizing configurations. This has 
strong consequences for the IDS at the bottom of the spectrum. An extreme case is given by 
the ID Bernoulli displacement model, i.e. /i = |^d„„^ + |^-d„a^, whose IDS satisfies the lower 
bound N{Eo + s) > C j In e. This singular behavior is the extreme opposite of a Lifshitz tail. 

(iii) The previous remark might mislead into expecting that the Bernoulli displacement 
model is not localized at low energy. Spectral localization for the ID random displacement 
model (at all energies and for arbitrary non-trivial distribution of the displacements) has been 
proven in [4] and [24] (using methods of [7]). These methods are completely different from 
what is available for ci > 1, and, in particular, do not require smallness of the IDS. One 
uses dynamical systems tools such as results on products of random matrices, in particular 
Furstenberg's theorem. In fact, one finds that the Lyapunov exponents are positive with the 
possible exception of a discrete set of energies. As far as dynamical localization is concerned, 
it might be violated at those critical energies. For examples of this see [16, 6]. Away from the 
critical energies, however, one also has dynamical localization. 

In the case of the one-dimensional Bernoulli displacement model, the energy Eq provides a 
new example of a critical energy. This is seen as follows: By the results of [2], for any e > 
almost surely there is a solution mq of H^u = Equ and C > such that ^ exp(— x~^/^~^) < 
|'Wo(3;)| < Cexp(x^/^~'"^). Using the lower bound and the standard reduction of order argument 
one shows that there is a second linearly independent solution ui which satisfies the same upper 
bound. This shows that the transfer matrices grow sub-exponentially. Thus, the Lyapunov 
exponent at Eq vanishes. 

(iv) An interesting open problem arises from cases of the random displacement model where, 
in the language used in the introduction, alternative (ii) holds, i.e. where EQ{a) vanishes 
identically. Examples for this (non-generic) situation can be constructed as follows: 

Let < ip E C°°(Ai(0)) be constant near the boundary (but not constant throughout 
Ai(0)). In the definition (1), (2) of the random displacement iJ^ model pick the single- 
site potential as g = Aip/ip. By construction, this leads to alternative (ii). Moreover, for 
every displacement configuration u a generalized ground state of H^ to -Eq = is given by 
= ^„ ip{- — n — u}n)Xhi(n) (here we think of ip as extended by a constant onto all of W^). 
Note that 1/C < < C for some C > uniformly in u. 

This leads to van Hove behavior of the IDS at Eq = 0, i.e. 

for E > 0, which follows with the same argument as provided for a closely related example in 
Section 3 of [20]. 

It would be interesting to know if this can generate non-trivial transport (and thus prevent 
dynamical localization). This is the case in dimension d = 1. Starting with 0, the reduction of 
order argument provides a second, linearly independent, solution which grows at most linearly. 
Thus, the transfer matrix also grows linearly. By Corollary 2.1 in [6], which is based on work 
in [8], this implies that the time-averaged moments \X\p of suitable solutions of the time- 
dependent Schrodinger equation are bounded below by CT^^~^^^'^. This rules out dynamical 
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localization in the sense of (55) for p > 5. Clearly, there are multiple obstacles to extending 
these methods to higher dimension. 

We finally remark that the methods of [8] and [6] do not suffice to obtain non-trivial trans- 
port under the sub-exponential growth bounds on the transfer matrix discussed in remark (iii) 
above. 
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